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Abstract. 

While many-particle entanglement can be found in natural solids and 
strongly interacting atomic and molecular gases, generating highly entangled 
states between weaJfly interax;ting particles in a controlled and scalable way 
presents a significant challenge. We describe here a one-step method to generate 
entanglement in a dilute gas of cold polar molecules. For molecules in optical 
traps separated by a few micrometers, we show that maximally entangled states 
can be created using the strong off-resonant pulses that are routinely used in 
molecular alignment experiments. We show that the resulting alignment-mediated 
entanglement can be detected by measuring laser-induced fluorescence with single- 
site resolution and that signatures of this molecular entanglement also appear 
in the microwave absorption spectra of the molecular ensemble. We analyze 
the robustness of these entangled molecular states with respect to intensity 
fluctuations of the trapping laser and discuss possible applications of the system 
for quantum information processing. 
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The concept of entanglement has evolved from being regarded as a perplexing 
and even undesirable consequence of quantum mechanics in the early studies by 
Schrodinger [1, 2] and Einstein [3], to being now widely considered as a fundamental 
technological resource that can be harnessed in order to perform tasks that exceed 
the capabilities of classical systems [4]. Besides its pioneering applications in 
secure communication protocols and quantum computing [5], entanglement has also 
been found to be an important unifying concept in the analysis of magnetism 
[6, 7, 8, 9], electron correlations [10] and quantum phase transitions [11, 12, 9]. Many 
properties and applications of entanglement have been demonstrated using a variety 
of physical systems including photons [13, 14, 15, 16, 17], trapped neutral atoms 
[18, 19, 20, 21, 22], trapped ions [23, 24, 25, 26, 27], and hybrid architectures [28, 29]. 
Entanglement has also been shown to persist in macroscopic [30, 31, 32, 33] and 
biological systems [34, 35]. Despite this significant progress, the theory of quantum 
entanglement and its technological implications are still far from being completely 
understood [4]. 

Trapped neutral atoms are regarded as a promising platform for applications 

of quantum entanglement due their relatively long coherence times [19], which can 
exceed those of solid state and trapped ion architectures by orders of magnitude 
[36]. Moreover, the sources of single-particle decoherence are well characterized in 
electromagnetic traps [19], and can be compensated using standard state transfer 
techniques [37]. In order to address individual atoms in an optical trap for coherent 
state manipulation, it is necessary to separate the particles from each other by a 
distance comparable to optical wavelengths [38, 39]. However, it is difficult to achieve 
entanglement between ground state atoms at such long distances, due to the short 
range nature of their mutual interaction. It is nevertheless is possible to enhance 
interactions between atoms in optical traps by either controlling the interatomic 
distance [40, 41, 42], or exciting atoms to an internal state that supports long- 
range interactions [43, 44, 45, 46, 47]. Using these methods, recent experiments 
have demonstrated the generation and characterization of entangled atomic states 
[18, 48, 21, 22], which are the first steps towards the study of many-particle 
entanglement and the development of quantum technologies using optically trapped 
particles. 

Quantum entanglement can also be studied using trapped polar molecules [49]. 
Arrays of polar molecules can be prepared in optical lattices with full control over 
the internal states including the hyperfine structure [50, 51, 52, 53]. Trapped 
molecules inherit the long coherence times of their atomic counterparts and the long- 
range dipole-dipole interaction between molecules offers a route for entanglement 
generation. Since the dipolc moment of freely rotating molecules averages to zero, 
proposals for molecular entanglement creation have involved the application of DC 
electric fields to spatially orient the dipoles [54]. One promising approach consists 
of placing the oriented dipoles in an ordered array using an optical lattice and 
performing entangling gate operations using microwave pulses, building on analogies 
with architectures for NMR quantum computation [55, 56, 57]. In order to overcome 
the complexity involved in controlling the "always-on" interaction between oriented 
dipoles, conditional transitions between weakly and strongly interacting states have 
also been proposed as a route to generation of intermolecular entanglement [58, 59, 60]. 
This approach has recently been demonstrated experimentally for cold atoms [21, 22]. 
Theoretical work has shown that entanglement can also be generated by coupling 
internal states with collective motional states in strongly interacting molecular arrays 
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[61, 62], analogously to methods developed for trapped ions [63]. In addition to these 
approaches for the controlled generation of pairwise entanglement between molecules, 
many-particle entanglement is also expected to emerge in the pseudo-spin dynamics 
of an ensemble of polar molecules with tunable interactions [64, 65, 66, 67, 68]. 

In contrast with previous approaches for generation of entanglement between 
dipolar molecules, the scheme proposed here does not involve the use of DC 
electric fields. Instead, we introduce here a method for deterministic generation of 
entanglement that uses strong optical laser pulses far-detuned from any vibronic 
transition. We consider closed-shell polar molecules in their ground rovibrational 
state, with each molecule individually confined in an optical trap in order to suppress 
collisional losses. We show that a single off-resonant laser pulse can mediate the 
entanglement of weakly interacting polar molecules separated by up to several 
micrometers. The degree of entanglement and the timescale of the entanglement 
operation are shown to have a well-defined dependence on experimental parameters 
such as the pulse intensity and duration. The laser parameters considered in this work 
are consistent with the technology developed to study molecular alignment in thermal 
gases [69, 70, 71, 72]. We note that entanglement of polar rigid rotors in strong laser 
fields has been considered before in the high-density regime [73, 74], where the dipole- 
dipole interaction energy is comparable to the rotational constant. The approach 
presented here allows for the generation of laser-mediated entanglement of rotors in 
dilute gases for the first time. 

The remainder of this paper is organized as follows. Section 1 reviews the 
rotational structure of closed-shell molecules in strong off-resonant optical fields. 
In Section 2 we analyze the generation of entanglement between two distant polar 
molecules due to the action of a single off-resonant laser pulse. The dependence of the 
degree of entanglement on experimental parameters is discussed in detail. In section 3 
we discuss two entanglement detection schemes, one based on Bell-type measurements 
for systems possessing single-molecule addressability and another scheme that employs 
microwave spectroscopy with only global addressing capability. In section 4 we 
investigate the effects of motional decoherence and show that entanglement in optical 
traps can be robust against this type of noise. We close with a summary and 
conclusions in Section 5. 

1. Molecules in far-detuned optical fields 

We consider closed-shell diatomic molecules in the vibrational and electronic ground 
state. The state of the molecules in the absence of external fields is represented by 
)iV, Mjv), which is an eigenstate of the rigid rotor Hamiltonian Hji ~ B^N"^ and 
Nz, where N is the rotational angular momentum operator and Nz its component 
along the space-fixed Z-asds. is the rotational constant. The interaction of a 
molecule with a monochromatic electromagnetic field E(r,t) = ^ [eii'(t)e*"* + c.c] 
whose frequency ui is far-detuned from any vibronic resonance can be described by the 
time-independent effective Hamiltonian [72] 



where Ep{r) is the space- fixed p-component of the positive- frequency field in the 
spherical basis and Q:p,p' is the molecular polarizability operator. For diatomic 




(1) 




Figure 1. Dimensionless rotational energy E/Bc of a molecule in the presence 
of a linearly-polarized CW far-detuned laser, as a function of the light-matter 
coupling strength Q.i = |i?opAa/4i?e: (a) Energies of the first six states with 
Mjv = (blue) and |M]v| = 1 (red). The states of the lowest doublet \g) = |0,0> 
and \e) = |1,0) define a two-level subspace. Be is the rotational constant, Ao 
is the polarizability anisotropy, and |i?oP = ^/2eoc, where I is the intensity of 
the laser. The notation |iV, Mjv) indicates that the rotational quantum number 
N is not conserved for f2j ^ 0. Mj^ is the projection of the rotational angular 
momentum along the laser polarization. 



molecules in a linearly polarized field, transforming the polarizability operator to the 
rotating body- fixed frame allows Eq. (1) to be rewritten as 

i^AC - |^(«|| + 2a^) + ^(aii - , (2) 

where v'^^l = {2,cos^ 6 ~ l)/2 is an element of the Wigner rotation matrix [75], Eq 
is the field amplitude for the selected polarization and 9 is the polar angle of the 
internuclear axis with respect to this. The polarizabilty tensor for diatomic molecules 
is parametrized by its parallel ay and perpendicular a±^ components, with a\\ > a±. 
The first term in Eq. (2) leads to a state-independent shift of the rotational levels 
and the second term induces coherences between rotational states |A^Mjv), according 
to the selection rules AA'^ = 0, ±2 and AAf^r = 0. Therefore the parity of rotational 
states in the presence of a far-detuned field is conserved. 

Ignoring the state-independent light shift (which contributes with just an overall 
phase to the eigenstates) and expressing the energy in units of B^, the single-molecule 
Hamiltonian H — Hji + Hj^c can then be written as 

H^N^-^n,vl^l{0), (3) 

where J7i = |i?oP(Q;|| — a±)/4:Bc is a dimensionless parameter that characterizes the 
strength of the light-matter interaction and is proportional to the field intensity 
Iq = ^ceol^'Op- In Fig- 1 wc plot the lowest eigenvalues of as a function of 
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Oj. The figure shows that for intense fields Cii » 10, the energy spectrum consists of 
closely spaced doublets, as first discussed in Rcf. [69]. The lowest doublet states \g) 
and |e) correlate adiabatically with the states \g) = |0, 0) and |e) = |1,0) in the limit 
fli — )• 0. Since the eigenstates of Hamiltonian in Eq. (3) have well-defined parity, the 
induced dipole moments (9|d|g) and (e|d|e) vanish, but the transition dipole moment 
(e|d|s() is finite for polar molecules, where d is the electric dipole operator. 

The light-matter interaction term Hac in Eq. (2) has been widely used to 
describe the alignment of polar and non-polar molecules in intense off-resonant fields 
[69, 76, 72]. From a classical point of view, the electric field of a strong off-resonant 
optical field polarizes the molecular charge distribution, inducing an instantaneous 
dipole moment. The field then exerts a torque on the rotating dipole that changes 
the angular momentum of the molecule, favouring the alignment of the dipole axis 
along the field polarization direction. However, the orientation of the dipole is not 
well-defined in AC electric fields. The degree of alignment for diatomic molecules is 
typically measured by the expectation value A = (cos^ 9) [72, 70, 71], with defined in 
Eq. (2). A is close to unity for aligned molecules. Adiabatic alignment in the presence 
of strong off-resonant laser pulses has been extensively studied both experimentally 
and theoretically [71, 72]. In adiabatic alignment experiments the laser pulse turn-on 
and turn-off times are long compared with the free rotational timescale t^ = h/Be- 
Under adiabatic conditions, the rotational motion of the molecules is described by the 
eigenstates of Eq. (3) with adiabatically varying values of Qi{t). 

In this work we consider molecules driven by strong off-resonant pulses that are 
adiabatic with respect to the rotational timcscales, but not necessarily adiabatic with 
respect to longer timescales such as the dipole-dipole interaction time between distant 
molecules (see below). 



2. Dynamical entanglement generation using strong laser pulses 

We now consider the dipole-dipole interaction between polar molecules in the presence 
of a strong off- resonant laser. The single- molecule Hamiltonian H = Hyi+Hag is given 
in Eq. (3) with intensity-dependent eigenvalues shown in Fig. 1. Using the two- level 
single-molecule subspace <Si = {15), |e)} the dipole-dipole interaction operator can be 
written as 

Kid = 7(l-3cos2e)[/dd(i?)x{|.9ie2)(ei,g2| + kiea) (.91.92! + H.c.} ,(4) 

where 7 = d~'^ {e\do\g)^ is a universal dimensionless parameter that depends on the 
external field strength and polarization, U^d = (P/R^ is the interaction energy scale, 
R is the intermolecular distance, Q is the polar angle of the intcrmolecular axis with 
respect to the laser polarization, do is the component of the electric dipole operator 
along the laser polarization and d is the permanent dipole moment of the molecule. 
At distances such that U^d/Bc ^ 1, the interaction operator Vdd does not mix the 
states \g) and |e) with higher field-dressed rotational states. 

The two-molecule Hamiltonian matrix H = Hi + H2 + Vdd in the subspace 
S2 = {1.91,92), |.9i, ^2), lei.92), |ei, 62)} can be written in two equivalent forms (up to 
a constant energy shift) as 

'H = ee (c\ci + cjc2 j + J12 {c\ + ci j (cl + C2) 

= |(4 + 4) + -/i2aif^i, (5) 
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where the operator cj = \ei){gi\ creates a rotational excitation on the i-th molecule, 
with the states jg,) and je^) equivalcntly represented by eigenstates of cr^ with 
eigenvalues — 1,+1, respectively, where ct^ (a = X,Y,Z) is a spin- 1/2 Pauli matrix. 
Ji2 = {eig2\Vdd\gie2) = (eie2|Vdd|5i52) is the exchange coupling energy, and Se is 
the splitting of the lowest doublet in Fig. 1. The eigenstates of H involving the 
single excitation sector are the symmetric and antisymmetric Bell states \'i!±) = 
2~^/^ {\g1e2) ± \e1g2)} with the eigenvalues E± = Se± J12. The ground and highest 
excited states can be written as 

!$_(«)) = cosa 15132) - sina |eie2) 
|$+(a)) = sina 15152) + cosa |eie2) ' 

with eigenvalues E± = ee± K, where K = i/e^ + J^^ . The states \^±{a)) are linear 
combinations of the remaining Bell states |$±) = 2-1/2 n^^^^^ ^ \e1e2)}. The mixing 
angle a is defined by tan(2Q;) = Ji^jsc- The states |<I'±(q;)) are separable in the limits 
a — )■ and a — > ±00. The ground state of the system is |$_(a)) for all values of a. 

Since the eigenstates of this two-molecule Hamiltonian are entangled for any finite 
value of the ratio Ji2/£o, we may consider the possibility of tuning the degree of 
entanglement by manipulating the transition energy with a strong off-resonant 
field. This corresponds to varying the effective magnetic field h = ec/2 for the spin 
chain Hamiltonian in Eq. (5). The possibility of preparing the states |$±(q:)) in Eq. 
(6) using strong continuous- wave (CW) off-resonant laser fields was first pointed out in 
Ref. [77]. However, since in practice the achievable intensity of CW lasers is limited, 
we consider here an alternative dynamical preparation of molecular entanglement using 
pulsed lasers. 

Polar molecules can be prepared in the rovibrational ground state |g) inside an 
optical trap [49]. A strong linearly polarized off-resonant field can then be used to 
bring the energy of the excited state |e) close to degeneracy with the ground state \g) 
by adiabatically following the energy level diagram in Fig. 1. In the presence of a laser 
pulse, both the dipolar coupling Ji2(t) and the excitation energy edt) become time- 
dependent. We take the initial two-molecule wavefunction as |*I'(0)) = 15152)- For this 
initial condition the state evolution is determined by the Hamiltonian sub-block 

V-i ° -^12(0 \ 



Ji2{t) 2e,it) 

with no participation of the single-excitation manifold since the Hamiltonian in Eq. 
(5) is block-diagonal. The state of the system is described by a superposition of the 
form 

|$(i))=a(t)|5i52)+6(i)|eie2). (8) 

Expressing the energy in units of the rotational constant B^. and time in units 
of = h/Bc, we can write the equations of motion m(r) = J{T)b{T) and 
ib{T) = J(T)a{T) + 2i?(T)&(r), which we integrate numerically using a standard 
Runge-Kutta-Felilberg method [78]. We have defined here the dimcnsionlcss energies 
J = Ji2/Bc, E = €,./ Be, and time r — t/t^. The dipole-dipole interaction timescale 
^dd = f>-/Udd depends on the intermolecular distance. The ratio between the rotational 
and interaction timescales tdd/^R is larger than unity for distances larger than the 
characteristic dipole radius (in atomic units) 

Ro = {dyB,)'^\ (9) 
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Figure 2. Evolution of the two-molecule concurrence C(p) under the action 
of a Gaussian off-resonant laser pulse with intensity profile = f^{t)Qot 

centered at t = 0. The intermolecular distance is -R = 10 Rq and the pulsewidth 
Tp = t(jd = lO^iR- Curves are labeled according to the value of the peak intensity 
Qq. The dashed line shows the envelope function of the pulse f{t). fdd is the 
dipole-dipole interaction time and = h/Bc is the rotational timescale. 



We solve the time-dependent Schrodinger equation by evaluating the energies E{t) 
and J{t) at each time step using an intensity parameter of the form fti{t) — [f{t)]'^no, 
for a Gaussian electric field envelope f{t) = e^'*/*") . We take to ^ to ensure 
adiabaticity with respect to the rotational motion. Under this condition we may 
extract E{t) from Fig. 1. The exchange energy Ji2{t) is evaluated using the 
instantaneous eigenstates \g{t)) and \e{t)) of the single- molecule Hamiltonian in Eq. 
(3). The parameter 7 varies in the range 1/3 < 7 < 1 as a function of increasing 
monotonically from its lower bound at fii = and reaching unity asymptotically as fii 
increases. The presence of a weak DC electric field in addition to the time-dependent 
laser field significantly changes this simple behaviour. We discuss the effect of a DC 
field in detail in Appendix A. In the following we shall consider the evolution of the 
system in the absence of DC electric fields. 

2.1. Tuning entanglement with a single laser pulse 

We consider here pulses that are non-adiabatic with respect to the interaction 
timescale t^d — (R/Ro)^ ^r- For a laser pulse that is adiabatic with respect to both 
^R and idd, an initial separable two-particle state would simply acquire a dynamical 
phase after the pulse is over and no net entanglement would be created in the system. 

We define the entanglement radius Rc as the intermolecular separation at which 
the dipole-dipole interaction energy U^d is equal to the energy of the transition 
\g1g2) \e1e2), i.e., 

i?c = (rfV2£e)'^'. (10) 
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Figure 3. Asymptotic two-molecule concurrence C{p) as a function of the 
intermolecular distance R (in units of Ro), long after the action of a Gaussian 
off-resonant laser pulse. For the distance R = 100 Ro, we choose the pulsewidth 
Tp = 10® tR (FWHM) and peak intensity Qq = 270, to obtain a maximally 
entangled state with C(p) = 1. tdd/^R = (R/Ro)^ is the dipole-dipole interaction 
time in units of the rotational timescale tji = h/Bg. 



For two molecules within this radius, mixing of the states \g1g2) and |eie2) is 
energetically allowed in the presence of a strong laser pulse. In the absence of DC 
electric fields, on account of the exponentially decreasing splitting of the doublet 
states as a function of the intensity parameter ili, the entanglement radius Re 
increases exponentially with Slj. For concreteness, the value f^i = 300 corresponds 
to Re ~ 3000 i?o, which corresponds to distances of several micrometers between 
molecules (see Appendix A). 

Let us consider a pair of polar molecules separated by a distance Rq ^ R < Re, 
where both molecules are initially in their rotational ground states, i.e., |5'(0)) = 
I (71(72) ■ The evolution of this system in the presence of a single Gaussian laser pulse 
is given by Eq. (8) and depends on three independent parameters: the intermolecular 
distance R, the pulse peak intensity Qq, and the pulsewidth Tp (FWHM). We use the 
binary concurrence C(p) = 2|a6| to quantify the degree of entanglement of the time 
evolved state \^{t)). The concurrence, which completely determines the degree of 
entanglement of pure binary states [4, 9], vanishes for separable states and is unity for 
maximally entangled states. Fig. 2 shows the evolution of concurrence for a pair of 
molecules separated by i? = 10 i?o under the action of a strong off-resonant Gaussian 
pulse. The pulsewidth Tp is chosen equal to the dipole-dipole interaction time tddj 
while the peak intensity fig is varied. Figure 2 shows that molecular entanglement is 
created in the presence of the laser pulse and reaches an asymptotic constant value 
when the pulse is over. We find that the qualitative behaviour of the system evolution 
is independent of R, and Tp, but that the actual value of the asymptotic concurrence 
depends strongly on the choice of these parameters. 

Fig. 3 shows how the asymptotic concurrence C(p) depends on the intermolecular 
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Figure 4. Asymptotic concurrence C{p) as a function of tlie peak intensity 
parameter ilj, long after the action of a Gaussian off-resonant laser pulse. The 
intermolecular distance is R = 100 -Rq- Data is shown for different pulsewidths 
(FWHM): Tp = 4dd (circles), Tp = 3tdd/4 (diamonds), Tp = fdd/2 (triangles), 
and Tp = tdd/4 (squares), tdd = 10®% is the dipole-dipole interaction time and 
t^ = h/Bg is the rotational timescale. 



distance R, or equivalently on the interaction time t^^, for fixed pulse parameters 
Tp — 10^ and Hq = 270. We have chosen the pulse parameters here to ensure 
that two molecules separated hy R = 100 i?o {R < Rc) become maximally entangled 
(C(p) = 1). For smaller distances R < 100 i?o, the asymptotic concurrence has an 
oscillatory dependence on R. For such distances the pulsewidth Tp is longer than 
the corresponding interaction time t^d- The system undergoes Rabi-typc oscillations 
between the states \g1g2) and |eie2) while the pulse is on. The oscillation stops when 
the pulse is over, giving the asymptotic concurrence shown in Fig. 3. For larger 
distances R > 100i?0j the concurrence decays monotonically with R, and eventually 
for R^ Re there is no entanglement. In this case the pulsewidth Tp is smaller than t^d, 
and the state population does not have time to undergo a Rabi cycle. Our calculations 
show that the behaviour of the asymptotic concurrence in Fig. 3 is independent of 
the choice of pulse parameters fio and Tp. The fast decay of the entanglement with 
distance is particularly useful for an array of molecules. By choosing the laser pulse 
parameters flo and Tp appropriately, it is possible to prepare highly entangled states 
between nearest neighbours only. 

The dependence of the asymptotic concurrence C{p) on the laser pulse peak 
intensity flo is shown in Fig. 4. Data are shown for a fixed distance R — 100 i?o 
and for different values of the pulsewidth Tp. For all values of Tp, the concurrence 
is negligibly small below an intensity threshold, here fla sa 70, whose value depends 
on the intermolecular distance R. Independently of the pulsewidth, the asymptotic 
concurrence increases with the intensity above this threshold until it reaches the 
maximum value (C(p) = 1). For a given distance R, the maximum concurrence is 
achieved at smaller peak intensities fio when the pulsewidth is equal to the dipole- 
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RbCs 


1.238 


441 


0.0290 


0.4 


6.4 


1.15 


KRb 


0.615 


360 


0.0386 


0.7 


3.7 


0.86 


LiCs 


5.529 


327 


0.1940 


3.8 


9.3 


0.17 


LiRb 


4.168 


280 


0.2220 


5.0 


7.3 


0.15 



Table 1. Molecular parameters for selected polar alkali-metal dimers: 7o is the 
laser intensity corresponding to Cli = {4:1t/c) loAay /2Bg = 1. Rq = {(P/B^)^/^ 
is the characteristic length of the dipole-dipole interaction and (r = h/Bg is the 
timescale of the rotational motion. Values of the polarizability anisotropy Aay , 
dipole moment d and rotational constant Be are taken from Ref. [79] . 



dipole interaction time tdd- After reaching the maximum value, the concurrence 
decreases with intensity as the population of the doubly excited state \e1e2) exceeds 
= 1/2 in Eq. (8). In the strong field limit flo 00, when R and Tp = t^d are 
held constant, the population is completely transferred from \g1g2) to \e1e2), with no 
net entanglement creation. 

The presence of an intensity threshold for the creation of molecular entanglement 
in Fig. 4 can be related to the notion of entanglement radius Re described earlier. 
For molecules within this radius, the mixing of the ground state \g1g2) with the two- 
excitation state \e1e2) is energetically favourable since the energy ratio Ji2/2£e = 
7(1 — 3cos^ 0)(i?e/-R)^ exceeds unity. When this energy ratio is less than unity, 
the state mixing is suppressed and the concurrence becomes negligible. For a given 
distance R and pulsewidth Tp, the intensity threshold thus occurs at values of Qq for 
which Re/R ~ 1. In Fig. 4, R^ w 100 iio for Qq = 130. 

2.2. Example: alkali-metal dimers in optical lattices 

Table 2.2 lists the laser intensity Iq of a traveling wave corresponding to a light- 
matter interaction parameter ili = 1 for selected polar alkali-metal dimers that have 
been optically trapped at ultracold temperatures [49, 53]. Predicted values for the 
polarizability anisotropy Aay and rotational constants for the rovibrational ground 
state are taken from Ref. [79]. For alkali- metal dimers, Iq is on the order of 10^ — 10^ 
W/cm^. This is well within the realm of feasibility, since continuous- wave laser beams 
with frequencies in the mid-infrared region (A ^ 1 fim) can have intensities on the 
order of 10* W/cm^ when focused to micrometer size regions [80, 81], while intensities 
higher than 10^° W/cm^ can be achieved using pulsed lasers. Strong laser pulses are 
routinely used in molecular alignment experiments, with pulse durations varying from 
less than a femtosecond to hundreds of nanoseconds [70, 72]. 

We now consider the interaction of pairs of polar molecules with a strong off- 
resonant pulse when the molecules are trapped in individual sites of an optical lattice. 
Typical experimental lattice site separations arc in the range Ul = 400 — 1000 nm 
[38, 82]. For most alkali-metal dimers in Table 2.2, these distances correspond to 
R ~ 10^i?o- The results in Figs. 3 and 4 therefore show that highly-entangled 
states of molecules in different lattice sites can be prepared using a single laser pulse. 
For example, two LiRb molecules separated hy ul = 730 nm can be prepared in 
a maximally entangled state by using a single Gaussian pulse with peak intensity 
/ = 1.35 X 10^^ W/cm^ and pulsewidth Tp = t^ = 150 ns. These laser parameters 
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can be achieved using current technology [70]. It is therefore possible to generate 

highly entangled states in currently available optical lattice realizations by choosing 
the appropriate combination of parameters flo and Tp, regardless of the molecular 
species. 

3. Detection of molecular entanglement in optical traps 

In this section we discuss how the alignment-mediated entanglement created between 
polar molecules in different sites of an optically trapped molecular array may be 
observed experimentally. We first show that the pairwise entanglement created in 
an ensemble of molecules as described in Sec. 2 gives rise to coherent oscillations in 
the microwave absorption line shape. Thus the global entanglement of the ensemble 
may already be detected by measurement of the linear spectral response as a function 
of frequency. We then outline how the time dependence of an initially entangled 
state generated by a strong laser pulse that subsequently evolves under the free 
rotational Hamiltonian may be tracked using correlations between local orientation 
measurements and a Bell inequality analysis [83, 84]. For pairwise entanglement 
of a pure state, this allows a direct measurement of the concurrence measure of 
entanglement for the initially entangled state. This second entanglement detection 
scheme; requires either single site addressing resolution in an optical lattice or 
individual trapping in separate dipole traps. Such addressability is now possible for 
trapped atoms [21, 22, 39] and is a subject of much experimental effort for trapped 
molecules. In contrast, the first approach is more amenable to current technology 
because it requires only global and not individual addressing. 

To show how these two detection schemes work, we shall consider explicitly an 
ensemble of molecules trapped in individual sites of a double-well optical lattice. Such 
lattices can be prepared by superimposing standing waves with different periodicity 
[85, 86, 87, 88, 89]. When the distance between two neighbouring double wells is a few 
times longer than the separation between the double-well minima, the alignment- 
mediated entanglement operation described in Sec. 2 can be designed such that 
only molecules within a single double- well become entangled. Separability between 
neighboring pairs is ensured by increasing the distance between adjacent double wells. 
We consider identical independent molecular pairs here for simplicity. In practice, 
inhomogeneities in the entanglement preparation step would lead to a distribution of 
concurrence values throughout the array. 

In the remainder of this section we discuss the detection of entangled pairs 
initially prepared at time f = by a strong laser pulse in the pure state |$o) = 
^015152) +bo\eie2) and show how we may measure the value of the initial concurrence, 
C(Po) = 2|ao6o|. For times t > 0, each molecule of the pair evolves under the free 
rotational Hamiltonian Hn (Section 1). The state component \e1e2) therefore acquires 
a relative dynamical phase which may modify time-dependent observables but does 
not change the concurrence. Our analysis will show that we can effectively extract 
the initial state concurrence C{po) from both the linear absorption spectrum and 
orientational Bell inequality measurements. 

3.1. Global entanglement measure in optical lattices 

It is well known that the macroscopic response of an ensemble of particles to 
an external field is affected by the presence of entanglement in the system [9]. 
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In particular, thermodynamic properties such the heat capacity and magnetic 
susceptibility have been established as entanglement witnesses for spin chains [9, 90]. 
In this section we will identify the signatures of entanglement on the AC dielectric 
susceptibility of a gas sample of J\f identical molecules. For simplicity we consider 
an ensemble of identical entangled pairs but the results can readily be generalized to 
many-particle entangled states. 

In the absence of DC or near resonant AC electric fields, an ensemble of rotating 
polar molecules is unpolarizcd. An applied electric field E(t) creates a polarization 
P{t). To lowest order in the field, this polarization is given by 

^ = y^ dt' {(d(i')d(t))o - (d(i)d(iO)o} • E(t'), (11) 

where (• • •)o denotes an expectation vahie with respect to the state of the ensemble 
in the absence of the external field. Typically the system is in a thermal state 
p = Z~^{l3)e-I^"\ where Hq is the field-free Hamiltonian, Z(^) = Trje"'^^"} is the 
partition function and = fc^T. For equilibrium states the autocorrelation function 
{A{t)B{t'))Q depends only on the time difference r = t—t'. As noted above, for analysis 
of the entanglement after the strong laser pulse is switched off, the Hamiltonian Hq is 
given by the two-molccule Hamiltonian T-L in Eq. (5) with Oj = 0. 

Given the polarization, Eq. (11), the microwave susceptibility for a thermal 
ensemble can be written as [91] 

X(u;) = -NPoifi) ^— ^, (12) 

where Po{f3) < 1 is the thermal population of the rotational ground state |0,0), and 
7e is decay rate of the rotational excited state |1, 0). The absorption spectrum is given 

by 

where A{lu) = Im{x(w)} and Fg = hj^, is the transition linewidth. 

Let us now consider the microwave susceptibility for an ensemble of entangled 
pairs initially prepared in the pure state |$o) = aol.91,92) +&o|'5i(;2)- Unlike the thermal 
case, the corresponding density matrix po — |<I'o)('I'o| describes a non-stationary state, 
with coherences that evolve according to Hq (in the absence of external perturbations) . 
In this case the response of the system to the field E(t) is given by Eq. (11) as for the 
thermal case, but the autocorrelation function (<l>o|d(i)d(i')|<l?o) now depends on the 
absolute values of the time arguments t and t', where these are defined with respect 
to a common initial time. 

The eigenstates of the coupled pairs in the limit Ji2/2eo ^ 1 are |$i) = \g1g2) 
with energy Ei = 0, I^-a) = 2-1/^ [15162) - |ei5f2)] with energy Ea = - J12, 
|*s) = 2-1/2 [|^^e2) + \e1g2)] with energy Es = Se + J12, and ($4) = [6162) with 
energy E4 = 2se (see Eq. (6)). The energetic ordering of the states [^'a) and |^'b) 
depends on the sign of J12. Using the non-stationary state i>o in the Kubo formula 
of Eq. (11), the microwave absorption spectra at frequencies u ~ losi = {Es — Ei)/h 
can be written as 

is 



A{uj) = M 



3h 



ml 



{UJSI - w)2 + 7I 
{u)si - w)2 -t-72_ 



+l«o&o| ,,., "\2 , .,2 > (14) 
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where A/p = Af/2 is the number of pairs, 7s is the decay rate of the state ^s- In 
the derivation of Eq. (14)wc have used the transition dipole moments (5's|d|$i) = 
V2{e\d\g) = ($4|d|*5), and (*A|d|$i) = = {<^i\d\^A}- The function jF^{t) 
contains the time dependence from the evolution of the entangled state under Hq 
and can be written as 

^c. (t) = e-T"* [(wsi - uj) sin (^41 (t) + 7s cos (j>4i (*)] , (15) 

where 041 (t) = uj4it — O^a is the free phase evolution of the two-molecule coherence, 
9ba is the relative phase of the two components of the initial state, defined by 
agfeo = |ao''o|e'^'"', and 741 is a decoherence rate introduced to account for dephasing 
channels. 

The amplitude of the time-dependent lineshape depends on the magnitude of the 
two- molecule coherence |ao6o| = C(po)/2. For a maximally entangled two- molecule 
state |$o) with relative phase 6ba = 0, the peak absorption (per molecule) at the 
resonance frequency w = wsi is 

~M " ¥ f cos(2a;esi)] • (16) 

The presence of dynamical peaks in the absorption or emission spectra is a 
general feature of wavepacket evolution that has been widely studied for single 
atoms and molecules [91]. More recently, the coherent oscillation of spectral peaks 
in the nonlinear optical response of molecular aggregates has been associated with 
entanglement between molecular units [35, 92]. Equation (14) shows that it is possible 
to identify entanglement in an ensemble of dipolar molecular pairs by measuring 
the linear absorption spectra. The procedure would be as follows. After preparing 
the system in an entangled state using a strong off-resonant laser pulse, a weak 
microwave field tuned near resonance with the lowest dipole-allowed transition would 
give an absorption spectrum whose line width shows damped oscillations at frequency 
a;4i = ABf./Ti. The presence of oscillations serves as an entanglement witness. 
Eq. (15) shows that the amplitude of this linewidth oscillation is proportional to 
the concurrence C{pq) = 2|ao6o| of the initially prepared state, while the decay of the 
oscillation depends on the decoherence rate 741. Measuring the amplitude of these 
oscillations can thus allow measurement of the pairwise entanglement between the 
dipolar molecules. 



3.2. Bell's inequality for orientation correlations 

Bell inequalities quantify the differences between quantum and classical correlations 

of measurements performed in different bases on quantum systems and provide critical 
tests of the incompatibility of quantum mechanics with local realism. Violation 
of a Bell inequality constitutes evidence of nonlocal quantum correlations such as 
entanglement between distant particles [93]. Not all entangled bipartite states violate 
the inequality, although all separable states do satisfy the inequality [94, 95]. For the 
case of entangled molecules in the presence of DC electric fields, it was recently shown 
that violations of Bell inequalities can be established [83, 84]. In the following we 
adapt and simplify the analysis in Ref. [83] to analize the orientational entanglement 
of polar molecules trapped in an optical double well lattice and prepared in the pure 
state |$o) = ao\gig2) + &o|'5i'52) by the action of a strong off- resonant laser pulse. 
We assume that the subsequent evolution is determined as in Sec. 3.1 by the field- free 
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rigid rotor Hamiltonian Hn, i-e., we neglect the small perturbation due to the trapping 
potential. 

The degree of orientation of a single molecule is given by the expectation value 
of the operator O = cos0 [71, 72], where is the polar angle of the internuclear 
axis with respect to the quantization axis. The orientation operator in the two- 
level basis Si = {\g) = |0, 0),|e) = |1,0)} can be written as O = ax/V^, with 
eigenvalues A± = ±l/\/3, corresponding to the molecule being oriented parallel (plus 
sign) or antiparallel (minus sign) to the direction of the quantization axis. For our 
proposed realization with molecules trapped in double well optical lattices, orientation 
measurements can be performed in a using laser-induced fluoresence [96] with single- 
site resolution. 

We consider the two-time orientation correlation function for a molecular pair 
E{ti,t2) = (Oi(ti) ® 02(^2)), where 6i{ti) = U} {ti)d{0)Ui{ti) [83, 84, 97]. The free 
evolution operator is given by U{t) = q-'^^r*/^^ where Hb. = B^az in the two-level 
basis. The orientation correlation vanishes for separable two-molecule states, but 
remains finite for entangled states. In particular for a pair of molecules initially in the 
state |$o) = o,o\gi92) + &o|eie2), the orientation correlation function is given by 

E{ti,t2) = ^C(po) cos (Wegil + i^egt2 + dba) , (17) 

where C{po) is the concurrence of the initial pure state po = |(f>o)($o|, 0ba the 
relative phase between the state components (see above) and the rotational frequency 
is uieg = 2Bg/h. The correlation function is invariant under particle exchange and 
symmetric around ^1=^2= 7''/2 for the relative phase 9ha = nn, with n an integer. 

Bell measurements can be divided into three steps [93] . First is the preparation of 
a pair of particles, typically spins, in a repeatable way. Second, an experimental setting 
is chosen independently for each particle. The setting for spins corresponds to the 
orientation of a Stern-Gerlach apparatus that measures the spin projections of particles 
A and B along the directions a and b, respectively. Finally, the correlation E{a, b) 
between the measurement outcomes for different sets of directions (a, b) are collected. 
For quantum correlation the Bell's inequality in the Clauser-Horne-Shimony-Holt form 
[98, 4] 

\E{d, b) + E{a, b') + E{a!, b) - E{a', b')\ < 2A^^ (18) 

is violated, where Amax is the maximum value of the measurement outcome. The 
quantum mechanical spin projection operator is a - a, with a = (o"x, o"y, ct^). For 
spin-1/2 particles Amax = 1- 

There is a one-to-one correspondence between Bell measurements based on spin 
orientations a and b and a scheme based on the free rotational evolution of molecules. 
In the two-state basis used here, the molecular orientation operator in the Heisenberg 
picture can be written as 0{Ta) = -t^e*'^^'^''/^ ax e"*'^^'^"/^ = a ■ a, where we have 

defined the orientation vector a = ( l/-\/3) (cos Ta, — sin Ta, 0), and Tq = 2Beta/h. The 
time evolution of the orientation operator 0{Ta) thus corresponds to a clockwise 
rotation of the orientation direction a from the positive X axis by an angle Tq in the 
XY plane. Therefore, choosing the time ta when to perform a molecular orientation 
measurement is equivalent to choosing the orientation of the Stern-Gerlach apparatus 
for the case of spin-1/2 particles. The two-time orientation correlator in Eq. (17) can 
thus be written as E{ta, tb) = {d-a(S>b-a), which is the form of the correlation function 
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Time (units of T„) 

Figure 5. Violation of Bell's inequality for molecular orientation correlations. 
The absolute value of Si{t) = E{0,0) + E(0,t) + E{t,0) - E{t,t) is plotted as a 
function of time for several states of the form |<I>) = |ci||(;i(;2) + |^'|e'®''" |eie2}. Each 
panel shows l^i | for three values of the concurrence: C = 1.0 (black line), C = 0.9 
(red line), and C = 0.8 (blue line). Panels (a) and (b) correspond to the relative 
phases 9^,a = and Oi,^ = 7r/4, respectively. E{t,t') is the two-time orientation 
correlation function. Time in is units of the rotational period Tji = irh/Ba- 



for spin systems. FoUowing the equivalence between spin orientation and rotational 
evolution, the magnitude of the quantity 

S = E{ta, 4) + E{ta,t[) + E{t'^, 4) - E{t'^, tl). (19) 

can then be used to test violations of Bell's inequality. For our purposes it is sufficient 
to set ta — tf, = and ij, = tj, = t in Eq. (19) and evaluate the absolute value 
of Si{t) ^ £(0,0) + E{0,t) + E{t,0) - E{t,t) using Eq. (17). In Fig. 5 we plot 
|S'i(t)| as a function of time for several parent states |$o) = ao|<?i52) + ^o|eie2) with 
different concurrences C{p) and relative phases Oha- The upper bound imposed by 
Bell's inequality over the |5i(t)| is 2A^„a^x = 2/3. For the states shown in Fig. 5, this 
limit is violated over a wide range of times within a rotational period Tr — irt^ . The 
violation of the classical bound serves as an entanglement witness. Most importantly, 
the figure clearly shows that the degree of violation of Bell's inequality depends on 
the concurrence C(po) of the entangled state. Therefore, once the signal is calibrated 
it should be possible to use the magnitude of Si{t) at a chosen time to quantify the 
molecular entanglement. 

We close with some comments on experimental feasibility of these measurements. 
The preparation of entangled pairs can be done using the methods described in Sec. 
2. An ensemble of identical pairs can be prepared to enhance the sensitivity of the 
correlation measurements. Performing orientation measurements in individual sites 
with laser-induced fluorescence [96] is significantly less destructive than femtosecond 
photodissociation measurements. Experimental violations of Bell's inequality have 
been established in a large number of experiments using photons [99, 13, 15, 14, 100], 
trapped atoms [101], superconducting junctions [102], quantum dots [103], and even 
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elementary particles [104], but to the best of our knowledge it has not been established 
with molecules. Our analysis shows that it is possible with current technology to look 
for violations of Bell's inequality for molecules in long-wavelength optical lattices or 
in separate dipole traps. 

4. Robustness of entanglement against motional decoherence 

Entanglement between distant molecules can be expected to decay in time due to 
relaxation and dephasing processes resulting from environmental perturbations. For 
entangled molecules in optical traps decoherence processes arise from their interaction 
with noisy external fields. Far-detuned optical traps, for example, are sensitive to 
laser intensity fluctuations and beam pointing noise, which can cause heating of the 
trapped atoms or molecules [105, 106]. Trap noise affects the precision of atomic 
clocks [107, 108] and also the dynamics of strongly-correlated cold atomic cnscnibk^s 
[109]. Additional sources of decoherence influence the dynamics of the system in 
the presence of static electric and magnetic fields [110]. In this Section we analyze 
the robiistncss of alignment-mediated entanglement of molecules trapped in optical 
lattices to fluctuations in the optical trapping laser fields. Our primary focus here 
is on motional decoherence in optical arrays, which is most sensitive to the effective 
lattice temperature. 

For an array of interacting polar molecules, the fluctuation of the dipole-dipole 
interaction energy Um{R) with the motion of the molecules in the trapping potential 
represents a source of decoherence for the collective rotational state dynamics. The 
vibrational motion of the molecules in an optical lattice potential can be represented by 
phonons interacting with the coherent rotational excitation transfer between molecules 
in different sites. Following Ref. [Ill] we write the Hamiltonian for a one-dimensional 
molecular array in the absence of static electric fields as 

T-L = eggclci + Jijclcj 

i id 

+ Yl ^^k&lak ^^J^lc, (afe + 4) , (20) 

k k 

where aj, creates a phonon in the fc-th normal mode with frequency ojk- The first 
and second terms determine the coherent state transfer between molecules in different 
sites, with site energy Cgg = 28^ and hopping amplitude Jij (evaluated at equilibrium 
distances). The third term describes the vibrational energy of the molecular center 
of mass in the trapping potential, which we assume harmonic as an approximation. 
In the absence of DC electric fields the phonon spectrum is dispersionless [111], i.e., 
ojk = ojQ. The last term represents the interaction between the internal and external 
molecular degrees of freedom, characterized by the energy scale 

\^j{coo) = -3Ji2 

where ujq is the trapping frequency of the optical lattice, is the lattice constant, 
^0 = \Jfil2m,ijjQ is the oscillator length, and f^^ is a mode-coupling function that 
satisfies the relation = —fji- 

We have omitted terms of the form (c|cj -|- H.c) in Eq. (20), since these only 
affect the dynamics of the system when Ju/eeg ^ 1- As discussed in Section 2.1, this 
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condition is satisfied only in the presence of a strong off-resonant pulse. However, the 

laser pulse width Tp is orders of magnitude shorter than the timcscalc of the oscillation 
of molecules in the lattice potential (rp <C w^^^). This separation of timescales allows 
us to neglect the coupling between internal and translational degrees of freedom under 
the action of a strong off- resonant laser pulse, even when J/sgg ^ 1. After the pulse is 
over, the coupling to phonons can become important when the timescale for internal 
state evolution h/Ji2 is comparable with I/loq. Under this condition the molecular 
array evolves according to the Hamiltonian in Eq. (20) over a timescale shorter than 
the molecular trapping lifetime rtrap ~ 1 s [53]. 

The Hamiltonian in Eq. (20) can be rewritten as ^ = Hs + ^-3 +^-33 using the 
unitary transformation cj^ = X^i'^iM^l- The Hamiltonian Tis — X^^^m^JjC^ describes 
the collective rotational states in terms of excitonic states = cj^\g) with energy e^^. 
The second term T-Lb = Tt-^^o X!fe '^I'^fc describes free lattice phonons, and the term 

^SB=Y,X';jlcA&k + &l), (22) 

describes the interaction of the excitonic system with the phonon environment. The 
interaction energy in the exciton basis is given by A^^ = J^ij '^ln'^jv^ij- The internal 
state evolution of the excitonic system depends strongly on the characteristics of the 
phonon environment. For low phonon frequencies < Jn/h the interaction energy 
A*' can become the largest energy scale in the Hamiltonian, and non-Markovian effects 
in the evolution of the system density matrix p{t) become important [112]. We assume 
here for simplicity that TnJo > J12, or more precisely {lo/ aLY{Ji2l'hL0a) < 1 [113] so 
that we are in a weak coupling regime. Note that cjq is determined by the trapping 
strength of the optical lattice and that both this and the dipolar interaction J12 can be 
tuned in this system to a far greater extent than is possible for Hamiltonians describing 
excitonic energy transfer in molecular aggregates [114]. In this weak coupling regime, 
the system evolution can then be described by a quantum master equation in the 
Born-Markov and secular approximations [112]| as p{t) = —{i/K) [Hg, p{t)\-\-V {p{t)). 

Let us consider the case of two interacting polar molecules coupled to a common 
phonon environment via the nonlocal term in Eq. (22). The dissipative dynamics of 
the system density matrix p{t) is determined by 

V{p{t)) = l,v[~^p{t)v[-^ - \lM^\p{t)}, (23) 

where = \^s){^s\ ± are projection superoperators, 70 is the pure- 

dephasing rate, and {A, B} denotes the anticommutator. The projection into the 
two-excitation eigenstate 'P2 = |eie2)(eie2| does not contribute in the absence of DC 
electric fields (see discussion in Appendix A). The single-excitation eigenstates are 
\^s) = 2-i/2(|ei32) + \g1e2)) and \^a) = 2-i/2(|ei,g2) - Igie^))- Equation (23) shows 
that for a system prepared in the pure state |$) = a|5i32) +b\eie2) we have ^{p) = 0. 
In other words, the two-molecule entangled states prepared using a strong laser pulse 
do not decohere due to the interaction with environmental phonons in the optical 
lattice, regardless of the strength of the coupling to the environment and the effective 
lattice temperature. This is a consequence of the nonlocal nature of the interaction 
with the phonon environment and implies that under these conditions, the states 
|^±) = [15152) ± [6162)] provide a basis for a decoherence-free subspace in which all 
pairwise entangled states may be defined. 

X Note that the secular approximation does not allow for coherence transfer [34] 
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We can understand the effects of motional decoherence on the entangled triparticle 

and many-particle states by estimating the full phonon decoherence rates, given by 
7/ii/,/t'i/'('^), with indexing the excitonic states. In Eq. (23) the pure dephasing 
rate is defined as 70 = 'yAA,AA{0) = jss,ss{0) = -7aa,ss(0) = -jss,aa{0). In 
the Born-Markov and secular approximations, dephasing and relaxation processes 
that lead to decoherence and entanglement decay occur at the rate 7/tj/,/t'i/'(w) = 
(l/n^) J^^dTe''^^{Bf,^{T)B^>„,{0)), where (i?^,.(r)i3^v'(0)) is the bath correlation 

function with B^,^ = J2k ^^i/(^fc + ^l)- Appendix B we use a classical stochastic 

model to approximate the bath correlation function under the infliience of random 
intensity fluctuations of the trapping laser. This procedure allows us to write the 
decoherence rates as 

7m-,m'-' i^) = ^ ["(w) + 1] [J^l,^'.' (w) - J^i,^',' (-w)] , (24) 

where n{ui) = (e'^''" — 1)^^ is the Bose distribution function and 

J^i,M^) = E^^-^^'-' (^) (^_^f)2 + ^2' (25) 

is the semiclassical spectral density for optical lattice phonons. In Appendix B we 

show that the broadening parameter can be written as l3 = kuJq, where the factor 
K > is proportional to the strength of the laser intensity noise. The trapping noise 
causes damping of the correlation function as {B^^{t)B^i,{Q)) oc e~''l*l cos(a;'t), where 
w' = -^/cjq — The bath autocorrelation time Tc is order The condition for the 
Markov approximation to hold is thus j3~^ <C h/Ji2- 

For fixed trapping parameters wq, and /3, this analysis shows that different 
molecular species can undergo very different open system dynamics, depending on the 
strength of the dipolar interaction between molecules in different sites. For instance, 
let us consider LiCs {d = 5.5 D) and KRb [d = 0.6 D) species as examples of 
molecules with high and low permanent dipole moments, respectively. For an optical 
lattice with = 1 /xm and noise-induced damping rate j3 = 100 Hz, the open system 
dynamics would have Markovian behaviour for KRb molecules {Ji_2/h = 10 Hz), but 
for LiCs molecules {Ji2/h = 1.4 kHz) the system dynamics can be expected to be 
non-Markovian. A very attractive feature of this trapped dipolar molecule array is 
that the transition between Markovian and non-Markovian dynamics can be studied 
experimentally for any molecular species by manipulating the laser intensity noise in 
order to tune the parameter /3 as in Ref. [115], or by changing the lattice spacing 
to manipulate J12. 

In the regime where the Markov and secular approximations are valid, we 
can estimate the phonon-induced decoherence rate j{cos) in Eq. (24) (with state 
indices removed for simplicity) at the characteristic system frequency cos = Ji2/^- 
For a lattice temperature such that TiLOg/k^^T <C 1 the decoherence rate scales as 
7(a;s) ~ 4.^^{Ji2/hf{lo/aLfH{us), with H{uj) = (w/wo)/3/[(a; - ujof + p% For 
experimentally realizable parameters 13 = 1 kHz, = 10 kHz and Ol = 500 nm, 
the decoherence rate for KRb molecules {bJs/ln = 0.13 kHz) is 7(^5) ^ 10~^ Hz, 
which is negligibly small compared with the typical loss rate of molecules from optical 
traps (7trap ~ 1 Hz) due to incoherent Raman scattering of lattice photons. We 
conclude that the entangled states of polar molecules containing double excitations 
can be robust to phonon-induced decoherence in optical lattice settings for which the 
weak coupling condition %ujq/Ji2 ^ 1 holds. 
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5. Conclusion 

In this work we present a scheme to generate entanglement in arrays of opticaUy 
trapped polar molecules. Starting from an array of molecules prepared in their 
rovibrational ground state, a single strong off-resonant laser pulse can be used to 
generate entanglement between molecules in different sites of the array. The strong 
laser field induces the alignment of molecules along its polarization direction during the 
pulse. For such laser alignment of polar molecules interacting via a dipole-dipole term, 
the energy ratio between the coupling and site energies ^12/^6 can be larger than unity, 
allowing generation of two-particle wavefunctions of the form |$) = a\gig2) + b\eie2) 
in the presence of the strong laser field. For \ab\ ^ 0, the laser alignment will thus 
induce entangled states, where the precise form of the resulting entangled state may be 
controlled by the duration and strength of the laser pulse. The subsequent evolution 
after the laser pulse is completed adds a dynamical phase to the entangled state but 
does not change the concurrence measure of the extent of entanglement. The proposed 
generation scheme does not depend on the number of coupled molecules and also holds 
for a many-particle system. Here for simplicity we have considered explicitly only the 
two-particle case. 

We emphasize that this alignment-mediated entanglement involving double 

excitation states is not possible with static electric fields. The rotational structure 
of an aligned molecule is such that the transition energy £e between the lowest two 
rotational states \g) and |e) becomes comparable in magnitude with the dipole-dipole 
interaction energy Jij = (gi.gj|Vdd|eiej), for molecules separated by distances of up 
to several micrometers. At such large distances the ratio Ji2/ee is negligibly small 
in the absence of DC electric fields and double-excitation transitions of the type 
15152) 16162) are energetically suppressed. 

We have demonstrated explicitly that the degree of entanglement in a molecular 
pair can be manipulated by tuning experimental parameters such as the laser pulse 
intensity and duration, as well as the intcrmolecular distance. We presented two 
methods to detect and measure entanglement in optical traps after the strong laser 
pulse is applied. The first approach requires only global microwave addressing of the 
molecular array. Here wc showed that the linear microwave response of an ensemble of 
entangled pairs contains a contribution to the absorption lincshape that is proportional 
to the amount of pairwise entanglement and that oscillates in time at a frequency of 
order B^jh^ where Be is the rotational constant. Measuring the absorption peak 
oscillations over this timescale would then allow the concurrence of the state to be 
determined. The second approach is based on measurements of molecular orientation 
correlations to establish violations of Bell's inequality. This method relies on the 
ability to optically address individual sites of a molecular array in order to perform 
laser-induced fluorescence measurements. Finally, we also analyzed the robustness 
of the strong field alignment-mediated molecular entanglement in optical arrays with 
respect to motional decoherence induced by fluctuations in the trapping lasers. 

The results presented in this work for a molecular pair can readily be generalized 
to larger molecular arrays, as indicated in the text of the paper. In this context, it is 
useful to recognize that the system Hamiltonian can be mapped into a quantum-Ising 
model with a tunable magnetic field, a model that has been widely used in the study of 
quantum phase transitions [9] . Furthermore, the form of Ising Hamiltonian describing 
the system is 2-local, which supports universal quantum computation when combined 
with the ability to implement arbitrary single-particle unitary transformations [116]. 
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Therefore, an array of optically-trapped polar molecules driven by strong off-resonant 

laser pulses provides both a test-bed for studies of quantum entanglement in many- 
body systems and a novel platform for the development of quantum technologies. 
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Appendix A. Molecules in combined off-resonant laser and DC electric 
fields 

In this appendix we describe the dipole-dipole interaction between polar molecules 
in combined presence of DC electric fields and strong off-resonant pulsed laser fields. 
We discuss how the addition of a DC electric field affects the entanglement creation 
scheme described in Section 2. 

Dipole-dipole interaction in combined fields 

Let us consider a polar molecule in its vibrational ground state, under the influence 
of a DC electric field and a CW far-detuned optical field. If the laser polarization 

is collinear with the direction of the DC electric field (space-fixed Z axis), the 
dimensionless molecular Hamiltonian H = + Hdq + H ac can be written in analogy 
with Eq. (3) as 

H = N^- xv^]l - In.vl^l (A.l) 

where A = dEz/B^ parametrizes the strength of the DC electric field. Ez is the 
magnitude of the DC electric field and d is the permanent dipole moment of the 

molecule. The rotational structure for Ez = and large laser intensities Vl\ consists 
of harmonically spaced tunneling doublets separated by an energy proportional to ili 
as shown in Fig. 1 of the main text. Each doublet is composed of states with opposite 
parity whose energy splitting decreases exponentially with fii. Due to this near 
degeneracy, a very weak DC electric field strongly couples the field-dressed doublet 
states, splitting their energy levels linearly with A [117]. The two lowest doublet 
states \g) and |e) for A ^ 1 correlate adiabatically with \g) w yfa\Q,{)) + \/6|l,0) and 
|e) « a/^IO, 0) — -x/allj 0) ^'^ ^ 0; with b and \NMn) is an eigenstate of H^. 

In the absence of DC electric fields the dipole-dipole interaction operator Vdd has 
only one non-zero matrix element Jij ~ {eigj\Vdd\giej) — (e^e^ |Vdd|.9i.9j) , defined in 
Eq. (4). In the presence of DC electric fields the parity of the rotational states is 
broken and the following matrix elements become finite: Vfj^ = {9igj\Vdd\9i9j), = 
(ejej|t>dd|ejej), and V^f = (e^g^ |Fdd|e4.gj) • The dipolar energies {Jij.V^f -.V^f .V^j'} 
determine the dynamics of interacting polar molecules in the regime where the energy 
Aceg for the transition \g) |e) is much larger than the dipole-dipole energy 
Udd — d'^/R^, where R is the intcrmolecular distance. In the regime Aeeg ~ Udd 
two additional dipole-dipole transitions become important: Aij = {eigj\Vdd\gigj) and 
Bij = {eigj\Vdd\eiej). These matrix elements couple the single excitation manifold 
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Figure Al. Dipole-dipole interaction energies J12, D12 = ~ '12 
B12 as a function of the intensity parameter Oi. Curves are labeled according to 
tiic DC electric field strength A = dEz/Bc- The DC and AC electric fields arc 
coUincar. Energy is in units of U^d = d? /B? and the intermolecular axis is taken 
perpendicular to the orientation of the fields. 



with the ground and doubly excited states, and vanish in the absence of DC electric 
fields. 

In analogy with the definition of Jij in Eq (4) we can write the dipole-dipole 
energies in units of ?7dd(l - cos^ 8) as V^f = V^j' = /ieMg, = /^e' = A^egA^g, 
and Bij = /ieg/ie, where /leg = d~^{e\dQ\g) is the dimensionless transition dipole, 
He = d~'^{e\do\e) is the dimensionless dipole moment of the excited state and 
l^g = d~^{g\dQ\g) is the dipole moment of the ground state. For the choice of rotational 
states used here we have /Zgg > 0, /i^ > and /ig < 0, which give Aij = —Bij > 0. It is 
convenient to define the differential dipolar shift Dij — V^^^ — Vf^^ = Hg{fie ~ Mg) < 
to describe the single-excitation dynamics [111]. We evaluate the dipole-dipole matrix 
elements using the eigenvectors of the single-molecule Hamiltonian in Eq. (A.l). In 
Fig. Al we show the dependence of the dipole-dipole energies Jij, Dij, Aij and Bij 
on the laser intensity parameter fij and the DC field strength parameter A. The 
figure shows that the exchange interaction energy Jij tends to zero at high intensities 
r^i 3> 10 in the presence of a perturbatively small DC electric field A <C 1. The energies 
Aij and Bij also vanish at high intensities. Only the diagonal dipolar shifts V^^j^ , Vij^ 
and V^j^ are finite in the high intensity regime for any non-zero DC field strength. 

Disadvantages for dynamical entanglement creation 

The presence of a DC electric field modifies the state evolution under the action of 
a strong off-resonant laser pulse in two ways. First, a static electric field strongly 
mixes the quasi-degenerate doublet states at high laser intensities (Fig. 1), resulting 




in a linear DC Stark shift that increases the energy sphtting Sc- The Stark spUtting 
significantly modifies the entanglement radius i?o — (rf^/2eo)^^'^, as shown in Fig. 
A2. The value of Rc increases exponentially with the laser intensity parameter ilj 
in the absence of DC electric fields, but has an upper bound in combined fields. 
The bound depends on the DC field strength A = dEz/B^, which determines the 
splitting of the states \g) and |e). For larger values of A, the intermolecular distance 
at which the dipole-dipole interaction between molecules becomes comparable with 
the Stark splitting becomes smaller. For the molecular species used in Table 2.2, 
A ~ 1 corresponds to Ez ^ 1 kV/cm. For such large field strengths, Rc ^ Rq ^ I 
nm for most alkali-metal dimers. Therefore, molecules in optical lattices with site 
separation R 10^ nm cannot be entangled using strong off-resonant fields when 
DC electric fields Ez 1 kV/cm are present. Figure A2 however shows that in the 
presence of stray fields Ez < 1 mV/cm (A < 10^^), alignment-mediated entanglement 
of alkali-metal dimers in optical lattices is still possible. 

Second, breaking the parity symmetry of the rotational states results in additional 
contributions to the dipole-dipole interaction such as already discussed. The matrix 
elements Aij and Bij mix the subspaces Si = {\g1e2), \e1g2)} and ^2 = {I5152), |eie2)}, 
the two-molecule state for the initial condition |$(0)) = 1.91^2) is given by |<l'(t)) = 
a{t)\gi92} + b{t)\eig2) + c{t)\gie2) + d{t)\eie2), with \ad\ ^ and \bc\ ^ 0. Therefore, 
for intermolecular distances R < R^ the two-molecule state evolution in combined DC 
and off-resonant fields no longer follows the simple two-state dynamics described in 
Section 2. 

Static electric fields also affect the dynamics of the entangled states after the 
laser pulse is over. Local system-environment coupling occurs in the presence of a 
static electric field [111]. The local interaction of a pair of molecules with the phonon 
environment is described by Hint ~ k{c\ci + C2C2){a + d^), with k cx I?i2- The 
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associated dissipator can be written as 

V'{p{t))=joV[-^^p{t)vi^^ - lloM^lpit)} 

+ 4y„V2p{t)V2-2jl,{V2,p(t)}, (A.2) 

where V2 = \e1e2) (6162! is a Lindblad generator that induces dephasing of the doubly 
excited state. Therefore the two-molecule entangled state |$) = a\gig2) + b\eie2) no 
longer belongs to a Decoherence-Free Subspace (DFS) with respect to the phonon 
environment, i.e. 'D'{p{t)) ^ 0. The dccohcrence rate 7q would depend on the 
magnitude of the dipolar shift Dij , which can be tuned by manipulating the strength 
of an applied static electric fic^ld and the intensity of the trapping laser. In addition 
to the phonon-induced fluctuations of the site energies in the presence of DC electric 
fields, the molecular energies also undergo fluctuations dnc to electric field noise, which 
acts as a global source of dccohcrence that can lead to entanglement decay as discussed 
for general bipartite and tripartite states in Refs. [110]. 



Appendix B. Model spectral density of optical lattice phonons 

In this appendix we derive the expression for the transition rate in Eq. 

(24) using a semiclassical model for the phonon environment in optical lattices. 
We start from the system-bath interaction operator in the exciton basis Hsb = 

S^i/ Sfe '^Jti/cjjCj, ^Ofc -|- d^^ and define the time correlation function C^i,^^>,y'{t) = 

(B^^(f)S^'^'(O)), where the bath operator B{t) in the interaction picture is given 

by4>^(t) = EfcA^, [akit) + ai{t) . 

The classical vibrational energy of the array can be written as H = (1/2) Qj.+ 
oj'^Q\, where Qfe = Yl^=i^jk\/T^Xj are the normal modes of vibration defined in 
terms of the displacements Xj from equilibrium and the molecular mass m. Promoting 

normal coordinates to quantum operators as Qk = i/^T^'^ (o,k + allows us to 

write the semiclassical bath operator B^^{t) = J2k -^^i^\/^^Qfi^)- '^^^ classicalhath 
correlation function can thus be written as 

Cciit) - E ^IX'-^' (^) {Qkmrncu (B.l) 

where we used the fact that different modes {k' ^ k) are uncorrelated. The classical 
bath correlation function is a real quantity, i.e., C*j(f) = Cci(t). 

The quantum bath correlation function (omitting system state indices) is defined 
as C(t) = {B{t)B{0)) and satisfies C*{t) = C{-t) [112]. The system transition rate 
is given by 7(w) = G{uj)/Ti^ where G{uj) = dre'"^'^ C (t) is a real positive quantity. 
Using the detailed balance condition G{—lo) = e~^^"G{io), where (3 = l/k^T, it is 
possible to write 

where G^(w) = dre"^'^Im{C(r)}. We use this expression to obtain a semiclassical 
approximation to the quantum rate 7(1x1). 

The approximation scheme consists on relating the antisymmetric function G^(u;) 
to the Fourier transform Gci(w) = e*'^^Cci(r) of the classical bath correlation 
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function in Eq. (B.l). Following Ref. [118], we use « {(ihuj / 2) G R{io) , and 

postulate the scmiclassical cIostuc Cu{t) — Cci(t).This procedure is known as the 
harmonic approximation. The approximate quantum transition rate is thus given by 

'>H = ;]2r^GciH. (B.3) 

The next step is specific to the system considered here. It involves the evaluation 
of the correlation function {Qk{t)Qk{0))ci from the classical equations of motion of a 
molecule in the optical lattice potential. For simplicity, we consider the potential to 
have the harmonic form V{x) = ^muifx'^, where uj^ is the frequency of the normal 
mode k. The most general form of the mode frequency is ujk = u}of(k), where 
Wo = {'2/h)y/VLEit is the trapping frequency as determined by the lattice depth 
Vl and the recoil energy Eji of the molecule. The function f{k) accounts for the 
dispersion of the phonon spectrum and is determined by the dipole-dipole interaction 
between ground state molecules in different lattice sites [111]. In this work we consider 
molecules in the absence of static electric fields, therefore the induced dipole moment 
vanishes and the phonon spectrum is dispersionless. For any k, the mode frequency 
ojk = ojQ thus depends on the trapping laser intensity II since Vl oc II [38, 49]. The 
laser intensity noise therefore modulates the phonon frequency luq and can lead to 
heating when the noise amplitude is large enough [105, 106]. The motion of a molecule 
in a fluctuating harmonic potential can be modeled by the equation of motion (for 
each k) 

Qk+tjl{t)Qk=0, (B.4) 

where = cjq [1 + a^(t)], and a^{t) is proportional to the relative intensity noise, i.e, 
am cx (I Lit) - (/o))/(/o). 

The equation of motion in Eq. (B.4) is a stochastic differential equation with 
multiplicative noise, for which no exact analytical solution exists [119]. Using a 
cumulant expansion approach, the equation of motion for the correlation function 
{Q{t)Q{0)) can be written as [119] 

^{Qmm + 2/3-|(g(<)Q(o)) + w;2(g(i)g(o)) = o, (b.5) 

where /3 = a^a;QC2/4 is an effective noise- induced damping coefficient and ojq = 
0Jq{1 — a^uiQCi) is an effective oscillator frequency which includes a noise- induced shift 
from the deterministic value loq. Equation (B.5) is valid for all times provided arc <C 1, 
where Tc is the noise autocorrelation time. The coefflcients ci and C2 are related to 
the noise autocorrelation function by 

poo 

ci= im^it - r)) sin(2wor)dT (B.6) 
Jo 

poo 

C2= / {m^{t-T))[l-cos{2u:oT)]dT. (B.7) 
Jo 

The effective damping constant can thus be written as ^ = {a'^u}^/8)[S{0) — S'(2a;o)], 
where S{uj) = J^^{S,{t)£,{t — t))c^'^'^'^ dr is the noise spectral density. The dependence 
of the damping coefficient on the spectral density at twice the natural frequency 
indicates that this is parametric dynamical process that can lead to heating {(3 < 0) 
when S{2ujq) > S{0). Here we assume that the static laser noise is dominant and use 
/3 > 0, which is satisfied for trapping lasers with approximate 1/f noise as in Ref. 
[105]. 
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The s olution to Eq. (B.5) is (Q(t)Q(0)> = (02(o))c-^l*l cos(a;'t), with uj' = 
a/wq — . Wc have assumed the osciUator is underdaniped (loq > and ignored 
the noise-induced frequency shift (wg = wo)- The mean square amphtude ((5^(0)) can 
be obtained by averaging over initial conditions using Boltzmann statistics. For an 
ensemble of identical one-dimensional harmonic oscillators we have ((3^(0)) k\jT /loq. 
Combining these results we can write the classical bath correlation function in Eq. 
(B.l) as 

By inserting the Fourier transform of Eq. (B.8) into Eq. (B.3) we obtain the 
semiclassical transition rate 

7p.,P v {^) = ^ [n{uj) + 1] (w) - J^l^^,,, (-w)] , (B.9) 

where n(uj) = (e'^'*'^ — 1)^^ is the Bose distribution function and we have defined the 
semiclassical phonon spectral density 

This approximate expression for J{uj) should be compared with exact phonon spectral 
density for an ensemble of free quantum oscillators J^v,n'v'{}^) = Ylik ^'liv^'li'v'^i'^ ~ 
uik), which also satisfies Eq. (B.9). 
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